Gluino production in ultrarelativistic heavy ion collisions and nuclear shadowing 
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In this article we investigate the influence of nuclear effects in the production of gluinos in nuclear 
collisions at the LHC, and estimate the transverse momentum dependence of the nuclear ratios 
= j^^/A-^^^ and Raa = I^^M'^^^p^- We demonstrate that depending on the 
magnitude of the nuclear effects, the production of gluinos could be enhanced, compared to proton- 
proton collisions. The study of these observables can be useful to determine the magnitude of the 
shadowing and antishadowing effects in the nuclear gluon distribution. Moreover, we test different 
SPS scenarios, corresponding to different soft SUSY breaking mechanisms, and find that the nuclear 
ratios are strongly dependent on that choice. 
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The main aim of the Large Hadron Collider (LHC), 
which is already running and soon will be in complete 
operation with 14 TeV, is to find the Higgs particle. That 
discovery may either confirm the Standard Model (SM) 
or open new windows towards new physics. Although the 
SM explain all experimental data except neutrino masses, 
there are many reasons to go beyond it. Some theoretical 
problems in the SM are: hierarchy problem, electroweak 
symmetry breaking, gauge coupling unification, etc T]. 
The Minimal Supersymmetric Standard Model (MSSM) 
is the simplest supersymmetric extension of the SM, be- 
ing a good candidate to Physics Beyond Standard Model 
[l,, In the MSSM, for each usual particle, one assigns 
a superpartner with oposite statistics: it means that for 
each boson there is a fermionic superpartner, and the re- 
verse in the case of fermions. In the strong sector, one 
has the so called supersymmetric QCD (sQCD), where 
besides the gluon (boson) and quarks (fermions), there 
are the corresponding superpartners: gluinos (fermions) 
and squarks (bosons) . On this model, the gluinos are the 
superpartners of gluons, they are color octet fermions 
and therefore they can not mix with other particles, as 
a result its mass is a parameter of soft SUSY breaking 
terms. Gluinos are Majorana fermions, expected to be 
one of the most massive MSSM sparticles, and therefore, 
their production is only feasible at very energetic ma- 
chines such as the LHC. The gluino and squark masses 
are still unknown parameters, but they cannot be smaller 
than around a half TeV, as predicted by several models 
for SUSY breaking. The "Snowmass Points and Slopes" 
(SPS) % are a set of benchmark points and parame- 
ter lines in the MSSM parameter space corresponding to 
different scenarios in the search for supersymmetry at 
present and future experiments (See [3] for a very nice 
review). The aim of this convention is reconstructing 



Scenario 


rrig (GeV) 


rriq {GeV) 


SPSla 


595.2 


539.9 


SPSlb 


916.1 


836.2 


SPS2 


784.4 


1533.6 


SPSS 


914.3 


818.3 


SPS4 


721.0 


732.2 


SPSS 


710.3 


643.9 


SPS6 


708.5 


641.3 


SPS7 


926.0 


861.3 


SPSS 


820.5 


1081.6 


SPS9 


1275.2 


1219.2 



TABLE I: The values of the masses of gluinos and squarks in 
the SPS scenarios. 



the fundamental supersymmetric theory, and its break- 
ing mechanism, from the experimental data. The dif- 
ferent scenarious correspond to three different kinds of 
models. The points SPS 1-6 are Minimal Supergravity 
(mSUGRA) model, SPS 7-8 are gauge-mediated symme- 
try breaking (GMSB) model, and SPS 9 are anomaly- 
mediated symmetry breaking (mAMSB) model (043)- 
Each set of parameters leads to different masses of the 
gluinos and squarks, which are the only relevant SUSY 
parameters in our study, and we shown their values in 
Tab. (HI . It will be shown below that the choice of SPS 
scenario affects the results for gluino production. 

Another aim of the LHC is to study the possible cre- 
ation and characterization of the so called quark gluon 
plasma (QGP), which is one of the predictions of the 
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FIG. 1: (Color online) Ratios Rj = xJa/Ax/n for the valence, sea quarks and gluons, predicted by the DS S], HKN [lol ]. 
EPS08 [ll[, and EPS09 [3 parametrizations at Q = 595 GeV and A = 208. The uncertainty band is shown for EPS09 nPDF's. 



Quantum Chromodynamics (QCD) (see e.g. Q). The 
heavy ion program at RHIC have brought many interest- 
ing results about the evidences of the QGP formation, 
which is in fact consistent with an almost perfect liquid 



0. Apart from the QGP, cold matter effects play also a 
very important role, changing the amount of interacting 
quarks and gluons in a given kinematic region. 



If the gluinos are found in proton-proton (pp) colli- 
sions {y/s = 14 TeV) at the LHC and if their masses are 
not much larger than 1 TeV, they might also be pro- 
duced in collisions involving nuclei - pA (proton-nucleus, 
y/s = 8.8 TeF) and AA (nucleus-nucleus, y/s = 5.5TeV) 
LHC modes. In this case, nuclear effects have to be 
considered in the searches for this supersimmetric par- 
ticles. One important initial state effect is the so called 
shadowing effect, which makes the parton distribution 
functions of the bound proton in a nucleus A (nPDFs) 
to be different from the usual PDFs in the free proton, 
ffix,Q'o) = Rfix,Ql)ffix,Ql), where Rf are the nu- 
clear modification ratios which parametrize the nuclear 
effects. There are several parametrizations of nuclear 
PDFs, based on different assumptions and techniques 
to perform a global fit of different sets of nuclear ex- 



perimental data using the DGLAP evolution equations: 
EKS98 1, DS i, HKN EPS08 [ul and EPS09 
(l^ . where the two later include different RHIC data 
for the first time. Also, EPS09 includes an uncertainty 
band around the central values. The typical x behav- 
ior of the nuclear modification ratios is the following: a 
supression for x < 10~^ (shadowing), followed by an in- 
creasement around 10~^ (antishadowing), again a supres- 
sion for X > 0.3 (EMC effect), and a bigger increasement 
when X approaches 1 (Fermi motion). The whole effect 
is usually called shadowing. 

To ilustrate how shadowing can influence the amount 
of partons in the nuclear medium, we show in Fig. [T]the 
results for a few nuclear modification ratios for the gluons 
(Rg), valence {uy, dy) and sea (ug, dg, s) quarks. Results 
for charm and bottom are not shown, since they are not 
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included in some of the parametrizations above (see sec. 
Ufor more details). The hard scale Q — 595 Gel^ is the 
gluino mass (SPSla scenario, shown in Tab. which is 
quite high. We did not include the EKS in our analysis, 
since this parametrization is not defined for such high Q 
values. Concerning Rg, the usual shadowing (suppres- 
sion) for very low x is present in all parametrizations, 
being very small for DS (5% suppression at x ~ 10~^, 
flat behavior), stronger for EPS08 (25% suppression at 
X ~ 10"'^) and moderate for HKN and EPS09 (15% sup- 
pression at a; ~ 10~^). However, for the processes con- 
sidered in this work, the small- a; region do not contribute 
(see below), and therefore we only show the relevant x 
domain. The shadowing is very smaller for x 10^'^, the 
DS and EPS being inside the EPS error band in many 
X regions (except for very high x). On the other hand, 
at larger x, antishadowing (enhancement) is present in 
EPS08, EPS09 {x < 10-1) and HKN (larger x), but not 
in DS. The behavior with increasing x is also different, 
being the growth steeper for EPS08, and smoothed out 
in later EPS09. Concerning the other parton species, 
Ruv ^ s-nd Rus ~ Rds for s-ll parameterizations 
except for HKN which show rather large differences. For 
moderate values of x, the HKN valence d and gluons 
are enhanced, the valence u is suppressed, while the sea 
HKN u has an enhancement followed by a suppression at 
larger x (EMC effect). There are many investigations on 
inclusive heavy quark, quarkonium and prompt photon 
production in central proton-nucleus and nucleus-nucleus 
collisions (See e.g. Refs. [13 21j), which try to help in 
constraining the nuclear PDFs from several observables. 
The variety of nuclear effects may also be relevant for 
gluino production, since there are contributing diagrams 
with both (anti)quarks and gluons in the initial state. 

In the case of gluino production, because of the large 
gluino masses, the values of x probed tend to be quite 
high (from x > 10"^ to almost 1), and then the antishad- 
owing, EMC effect and even Fermi motion may be im- 
portant (depending on the kinematic region and nuclear 
PDF), which may enhance the gluino production rate 
compared to that obtained from single nucleon collisions 
at the same energy. Therefore, whereas the smaller cen- 
ter of mass energy (5.5 TeV (AA) and 8.8 TeV (pA)) will 
reduce the gluino production rates (compared to 14 TeV 
(pp)), there may be an enhancement due to the amount 
of quarks and gluons on the nuclear medium compared 
with the nucleon parton distributions on a single proton, 
due to high density nuclear effects. In this work we inves- 
tigate whether this enhancement /suppression is present 
or not. 

This article is organized as follows. The basic formulae 
to calculate gluino production are presented in section H] 
Our results for gluino produced in nuclear collisions at 
the LHC are presented in section [TTl followed the the 
conclusions. 




(c) 

FIG. 2: Lowest order Feynman diagrams for gluino produc- 
tion: (a) quark-antiquark initial states and (b) gluon-gluon 
initial states (double gluino production); (c) quark-gluon ini- 
tial states (squark-gluino production). 



I. GLUINO PRODUCTION IN PP COLISIONS 



In order to make a consistent comparison and for sake 
of simplicity, we restrict ourselves to leading-order (LO) 
accuracy, where the partonic cross-sections for the pro- 
duction of squarks and gluinos in hadron collisions were 
calculated at the Born level already quite some time ago 
(2^ . The corresponding NLO calculation has already 
been done for the MSSM case [11], and the impact of the 
higher order terms is mainly on the normalization of the 
cross section [2^, which cancels out in the ratios. 

The contributing LO diagrams for inclusive gluino pro- 
duction in proton-proton collisions are gg ^ gg, qq ^ gq 
and the Compton process gq — > gq, (shown in Fig. 
where one has to be carefull in including Feynman rules 
for Majorana particles [23 ). 

Incoming quarks (including incoming b quarks) are as- 
sumed to be massless, such that we have Uf = 5 light 
flavours. We only consider final state squarks correspond- 
ing to the light quark flavours. All squark masses are 
taken equal to rriq (L-squarks and i?-squarks are there- 
fore mass-degenerate and experimentally indistinguish- 
able.) We do not consider in detail top squark production 
where these assumptions do not hold and which require 
a more dedicated treatment [1^. 

The invariant cross section for single gluino production 



can be written as 



da 



( C+cosf \ 
\^ 2 sine ) 



dt 



(ij -> gd), (1) 



where fij are the parton distributions of the incoming 
protons and ^ is the LO partonic cross section '22] for 
the subprocesses involved. The identified gluino is pro- 
duced at center-of-mass angle 9 and transverse momen- 
tum pt, and x± = The Mandelstam variables of 
the partonic reactions ij — >■ gg,gq are then 
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where m-g and are the masses of the final-state partons 
produced. The center-of-mass angle 9 and the differential 
cross section above can be easilly written in terms of the 
pseudorapidity variable r\ = — lntan(0/2), which is one 
of the experimental observables. 

Predictions for gluino production in pp collisions at the 
LHC (-^/i = 14 TeV), in all SPS scenarious, are shown in 
a former work , where there is a huge difference in the 
magnitude of px distributions for different SPS points, 
making it possible to distinguish between some different 
SUSY breaking scenarios. We can ask if the same occurs 
in nuclear processes, and answering this question is also 
a goal on this article. 



II. GLUINO PRODUCTION IN NUCLEAR 
COLLISIONS 

Let us now focus on gluino production in nuclear col- 
lisions. The calculation is done as explained in the pre- 
vious section, replacing the parton distributions in the 
free nucleon (ff in Eq. (H])) by the corresponding nu- 
clear parton distributions (for the proton PDF we 
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FIG. 3: (Color online) Transverse momentum dependence of 
the nuclear modification factor RpA for inclusive gluino pro- 
duction in pA collisions at the LHC {y^ = 8.8 TeV, jjyj < 2.5), 
for distict nPDFs. 



use the CTEQ6L1 27]). The nuclear effects are then 
studied by comparing the different nPDF's available (for 
consistency, we use the LO version of all nPDF's). To be 
sure that the nPDF's are within the regions of validity, 
we have used Q = nig as the hard scale (as done in [2j). 
Another possible choice, a px running Q = nig +pt scale, 
would push some of the nPDF's outside the region of vali- 
didy (EPS08 and EPS09 are frozen in Q = 1000 GeV for 
values above that scale, whereas DS is not valid in that 
region). For this reason, the DS could not be considered 
in the SPS9 scenario (see Table I) , with extra large gluino 
masses. To start with, we consider the SPSla scenario 
as the first (most optimistic) choice of gluino and squark 
masses. 

In Fig. [3] we show our results for the transverse mo- 
mentum dependence of the nuclear modification factor 
defined by 



R 



■pA 



d'^a{pA)_^^fa{pp) 



(4) 



drjdpT drjdpT 

for gluino production in proton-nucleus collisions at the 
LHC (Vs = 8.8 TeV). For lower pr, the DS and EPS08 
nPDF's are inside the EPS09 uncertainty band, with al- 
most no nuclear effect, RpA ~ 1- For pr > 500 GeT^, the 
EPS's starts to be slightely suppressed (increasing with 
Pt), whereas the DS starts to be slightely enhanced (in- 
creasing with Pt)- For the HKN distribution, there is a 
larger enhancement of 10%, increasing slowly with pT- 
This means that the correct amount of (anti)shadowing 
is undefined. In fact, as pT grows the probed values of 
X increase, and the EPS nPDF's enter the EMC region, 
whereas this effect does not appear for the other nPDF's. 

In Fig. 21 we present a similar analysis for the trans- 
verse momentum dependence of the nuclear modification 
factor defined by 



Raa 



d^a{AA)_^^2d^a{pp) 



drjdpT 



drjdpT 



(5) 



for gluino production in nucleus-nucleus collisions at the 
LHC {y/s = 5.5 TeV). In this case, the nuclear effects are 
amplified because of the presence of two nuclei. Besides, 
the probed values of x are pushed into very high- a; due 
to the smaller center of mass energy. Indeed, the EPS 
supression increases with px in a stronger way than in 
the pA case (around 15% for higher pr). The DS nPDF 
has an enhancement pattern, increasing with px, which 
shows that this distribution has reached the Fermi mo- 
tion effect in the very right side of Fig. 1. The enhance- 
ment is also larger for the HKN, above 20% with a very 
tiny increasement with px- It seems that, if the latest 
EPS09 nPDF is the more correct distribution, the gluino 
production will be slightely suppressed compared with pp 
collisions at the same energy, whereas the DS and HKN 
suggests that there will be some enhanced production of 
gluinos in nuclear collisions. 

The possible increasement of the gluino production 
rate in nuclear collisions (compared with pp collisions 
at same energy) shown above is in fact too low to re- 
ally improve the small feasibility of detecting the gluinos 
when going from pp to pA and A A. In fact, the higher 
hadronic activity in nuclear collisions make the detec- 
tion of gluinos more difficult, and the smaller CM energy 
available produces a smaller number of gluinos compared 
to 14 TeV pp colisions. The expected luminosity to be 



reached in the AA collisions (£ 



[28| is seven orders of magnitude smaller than in the pp 
mode {Cpp « 10'^^ cm~^s~^), and this is the main limita- 
tion to detecting nuclear gluinos (they will be produced 
but will hardly be seen). In the pA mode, one expects 
a luminosity of CpA ~ 7.4 x 10^^ cm~^s~^ [1^, which 
becomes 7Apb~^ assuming a full LHC year 10^ s (one 
usually considers a month ion running time 10^ s) in the 
ion mode. With only our LO estimation, and consider- 
ing the more suppressed EPS09, one would than obtain 
around 31 gluinos produced in the pA mode for the pr 
integrated region, so statistics is really limited. It has 
been suggested that the pA luminosity could eventually 
be upgraded to CpA ~ 10^^ cm~^s~^ [30|, in this case our 
estimate would increase to 430 gluinos in one year run. 
For more realistic estimates, the NLO correction would 
still increase the cross-sections for the various production 
processes by up to a factor of less than two [23] . 

Not only the nuclear shadowing but also the SUSY 
breaking parameters affect the nuclear ratios. This de- 
pendence is indirect, since the gluino and squark masses 
(TOg, rriq), are the only parameters that really affect the 
results, but these masses are consequences of the differ- 
ent SUSY breaking parameters in the different SPS sce- 
narious. This is shown in Fig. [5l where different SPS 
scenarious give different absolute values for the Raa nu- 
clear ratios (this can be seen by comparing for example 
the starting point of each curve) . The pr growth for the 
DS nPDF is even more steeper for the higher mass SPS 
scenarious (higher x). For the SPS9 scenario, the results 
should not be trusted, since most parametrizations are 
not valid in that region: the HKN predicts an enhance- 
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FIG. 4: (Color online) Transverse momentum dependence 
of the nuclear modification factor Raa for inclusive gluino 
production in AA collisions at the LHC (-^s — 5.5 TeV, 
i 77 1 < 2.5), for distict nPDFs. 



ment essentially constant with px, and the frozen EPS's 
supression decreases with px- Because of the odd in- 
terplay of nuclear effects and SUSY breaking scenarios, 
one needs to put better constraints on the nuclear PDF's 
before describing precisely gluino production in nuclear 
collisions. Conversely, the discovery and measurement 
of the gluino and squark masses will be important in 
the searches for sparticles produced in nuclear collisions, 
taking into account the correct nuclear effects which also 
depend on the those masses. 



III. CONCLUSIONS 

To conclude, in this work we studied the nuclear effects 
in pA and A A gluino production at the LHC. We have 
shown different results of enhancement or suppression de- 
pending on the nuclear PDF, the effects being smaller in 
pA interactions and larger in nuclei collisions. Gluinos 
will probably be copiously produced in the pp channel. 
Once the details of gluino production are known in pp 
interactions, studying this final state in pA and AA colli- 
sions could give unprecedented constraints on the nPDFs 
in a heretofore unexplored region of Q^. One could use 
the higher energy to get a good measurement of gluino 
production and search for deviations from that in the 
measurable px range for pA and AA to measure quark 
and gluon shadowing at very high scales where nothing 
at all is known about it. Uncertainties on the nPDF's 
(and cold matter effects in general), and on the SUSY 
breaking scenarios (which give different masses for the 
gluinos and squarks) has to be disentangled in the future 
searches. For heavy nuclei collisions, where its expected 
the formation of the quark gluon plasma, it may appear 
other channels where gluino is produced. Here we only 
investigated cold matter effects, namely the shadowing 
of the nuclear distributions. If gluinos are discovered in 
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FIG. 5; (Color online) Transverse momentum dependence of the ratio Raa in single gluino production at the LHC {^/s — 
5.5 TeV), for different choices of nuclear parton distributions: DS HKN O EPS08 and EPS09 [Q, in different SPS 
scenarios. 



pp collisions at LHC, they will also be there for pA and on a further understanding of the correct nuclear effects. 
AA. However, the ability to search for them will depend 
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